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Q , Abstract 

Type OA string theory in the (2, 4k) super confer mal minimal model backgrounds and 
, the bosonic string in the (2, 2k — 1) conformal minimal models, while perturbatively 

Qj I identical in some regimes, may be distinguished non-perturbatively using double scaled 

• matrix models. The resolvent of an associated Schrodinger operator plays three very 

important interconnected roles, which we explore perturbatively and non-perturbatively. 
On one hand, it acts as a source for placing D-branes and fluxes into the background, while 



\^ • on the other, it acts as a probe of the background, its first integral yielding the effective 



force on a scaled eigenvalue. We study this probe at disc, torus and annulus order in 
perturbation theory, in order to characterize the effects of D-branes and fluxes on the 
matrix eigenvalues. On a third hand, the integrated resolvent forms a representation of 
a twisted boson in an associated conformal field theory. The entire content of the closed 
string theory can be expressed in terms of Virasoro constraints on the partition function, 
which is realized as wavefunction in a coherent state of the boson. Remarkably, the 
D-brane or flux background is simply prepared by acting with a vertex operator of the 
twisted boson. This generates a number of sharp examples of open-closed duality, both 
old and new. We discuss whether the twisted boson conformal field theory can usefully 
be thought of as another holographic dual of the non-critical string theory. 



-'^Also: Visiting Professor at the Centre for Particle Theory, Department of Mathematical Sciences, University 
of Durham, Durham DHl 3LE, England. 
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1 Introduction and Discussion 



In recent years, it has become clear that the separation between open and closed string de- 
scriptions of certain physical phenomena is very much an artifact of perturbation theory. This 
means more than just strong/weak coupling dualities linking different string theories, but also 
includes transitions within a given theory from an open string description {e.g., involving D- 
branes) to a closed string one {e.g., involving only fluxes), as the strength of a coupling changes. 
Gauge/string dualities are the best known examples of this type of phenomenon, and for that 
reason alone it is of interest to learn more about how to describe them. For this, it is of course 
necessary to go beyond perturbation theory, and furthermore, to develop a far richer picture 
than that which is accessible using only simple strong/weak coupling duality. We need to be 
able to define the theories in question at arbitrary values of the couplings in order to address 
the many subtle issues which can arise. 

Matrix models in the double scaling limit [1-3] defined for us the earliest examples of exactly 
solvable string theories with readily accessible non-perturbative information. These types of 
models have recently been interpreted as examples of open/closed dualities in a very direct 
sense [4] , but also describe more subtle of interesting open and closed string phenomena, as we 
shall recall. In the direct sense, they represent certain decoupling limits of the world-volume 
theories of D-branes in non-critical string theories [5-7] , and at large — as in AdS / CFT — 
they reconstruct their parent string theories in specific spacetime backgrounds. These models 
are particularly interesting in that they are not supersymmetric, the D-branes are generically 
unstable, and as such their open string modes are tachyonic. The matrix models in fact supply 
us with perfect models of tachyon condensation [8,9], allowing for the unified description of 
D-branes, closed strings, and the transitions between them [10-15]. For example, the picture 
of D-brane decay as the process of rolling down a tachyon potential to find a new minimum 
representing closed strings always suffered somewhat from at least one conceptual discomfort: 
Once one arrives at the bottom of the potential, one still has the open string variables: Where 
are the closed strings? Without the description of the closed string theory to which one is 
supposed to have arrived, the picture drawn in figure ^a) is really only a metaphor^. In 
the double scaled matrix model approach, such diagrams are refined into more precise ones 
such as that depicted in figure ). The matrix model potential (after double scaling) is the 
tachyon potential. In the eigenvalue description, at large A^ the eigenvalues form a droplet 
(or collection of droplets) in one (or more) of the minima of this potential. The collective 
-^See e.g. ref. [16] for a discussfon. 
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Figure 1: The picture of tacliyon condensation: (a) Incomplete in general, due to a lack of a description 
of the endstate as closed strings, and (b) Complete in the matrix model description, since the sea of 
eigenvalues at the minimum is a holographic description of the closed string background. 

excitations of the eigenvalue constituents in these droplets are closed string physics, while the 
process of moving an eigenvalue out of the collective and up to a extremum of the potential 
represents the creation of a D-brane, whose tension is given by the energy difference between 
the sea of collective eigenvalues and the extremum. 

Furthermore, it was shown some time ago how to describe, with matrix models, string vacua 
containing a specific number (say, F) of background D-branes, both for non-perturbatively 
unstable models [17,18] and non-perturbatively stable ones [19,20]. We will be examining these 
definitions more in this paper. In each set of models, the region of large positive cosmological 
constant, /i, in each model give identical perturbation theory, and can be identified as having F 
background D-branes. The non-perturbatively stable models allow one to continue to a second 
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perturbative region, that of large negative whose interpretation was unclear at the time the 
models were first presented. Recent work [21] has shown that these non-perturbatively stable 
models are actually descriptions of the type OA string theories, and that the large negative /i 
perturbative regime has an interpretation in terms of there being F units of R-R flux. This is 
an example of an open/closed transition happening as a coupling (//) in the theory changes. 
(See refs. [22,23] for generalizations to the rest of the minimal models for type OA.) 

In this paper we wish to examine in detail some of the key components of these matrix model 
definitions, and sharpen the understanding of many of the features of the models. A key role is 
played by the diagonal of the resolvent R{z, As) of the Schrodinger operator H — —d^ + u{z), 
(As is a spectral parameter) which enters in three logically distinct (but interconnected) ways: 
(a)li acts as a source for placing D-branes and fluxes into the background, as can be seen from 
the way it enters the matrix model and the resulting string equations, (h) It acts as a probe 
of the background, since its first 2;-integral (after a Laplace transform) is the macroscopic loop 
expectation [24] operator. Its probe character also follows from the fact that its first 2;-integral 
yields the effective force on a scaled eigenvalue [25] . In order to understand better the physics 
in the various regimes, we switch on the background, corresponding to including the resolvent 
in the string equations with strength F, and then we work in each perturbative regime, using 
the integrated resolvent as a probe to see the effects of the background. We study this probe 
at disc order in large positive // perturbation theory, and to annulus and torus order in large 
negative // perturbation theory, in order to characterize the effects of D-branes and fluxes on 
the matrix eigenvalues. 

(c) The integrated resolvent forms a representation of a twisted boson ^p{z) in an associated 
conformal field theory. The entire content of the closed string theory can be expressed in terms 
of Virasoro constraints [26, 27] — the microscopic loop equations — on the partition function 
which is realized as wavefunction in a coherent state of the boson. The creation modes of the 
boson are closed string operator insertions in the string theory. The Virasoro constraints in the 
presence of F were worked out in ref. [19] and refined in ref. [20].^ It was observed in ref. [20] 
that the insertion of the D-brane (or flux) background in this language corresponds to just 
acting with a vertex operator'^ of the twisted boson, Vr =: exp(— F93((t)/2) : The parameter a 
corresponds to the boundary cosmological constant. 

This latter observation is particularly interesting since it makes extremely explicit (if it were 

not manifestly so already) the fact that in a fully non-perturbative framework, such as this 

^See also ref. [28], but note the comments on it in the discussion section of ref. [20]. 
•^It appears that this fact was rediscovered many years later [29]. 
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one, structures that we interpret as open strings in some perturbative regime can be introduced 
in entirely closed string terms, built out of a specific combination (derived in refs. [20,28]) of 
closed string operators'* . Here we see it succinctly in terms of the vertex operator. This is 
of course one of the biggest lessons to be taken away from all of this: Open/closed transitions 
and strong/weak coupling dualities exist because the physics can be specified independently 
of the perturbation theory (this is also the central lesson of M-theory). When we can get a 
handle on the right variables whose definition is not rooted in the world-sheet expansions, we 
can describe the physics in a number of ways, involving open or closed strings. 

This also leads one to wonder whether the twisted boson framework should be taken as an 
equivalent starting point for discussion of the non-critical string theory in various backgrounds. 
Is there a sense in which it is another (holographic) dual of the theory, or does it rise above 
that, being the parent theory within which holography and duality become manifest? This 
would go beyond the recent interesting statements in the literature [31-33] that the complex 
curve described by the analytically continued effective potential is simply an alternative target 
space. Instead, it would suggest that the twisted boson should be given dynamical meaning^. 
Specific solutions for the wavefunctions correspond to certain non-critical string backgrounds 
we already know, and so we should expect that the dynamics of can give us more insight into 
non-perturbative dynamics of string vacua, including perhaps time dependence, if we can find 
analogous structures in higher dimensional string theories. 

There is a fair amount of review material in this paper for which we make no apology. It sets 

up the language and conventions, and reminds the reader of older (and often highly relevant) 

literature. The following outline will help the well-informed reader avoid the parts they wish 

to: In section |21 we recall the string equations for the bosonic and type OA systems, noting 

that the former are perturbatively equivalent to the latter in the large positive cosmological 

constant regime. Section IHl is a detailed review of the double scaled matrix model origins of 

the bosonic string equations, while section |3] reviews the details of how to extract the effective 

(tachyon) potential for eigenvalues in the scaling limit. The non-perturbative effects of D- 

branes (instantons) are recalled. Their tensions are computed from the effective potentials for 

^This observation that one can move back and forth between closed and open string minimal model back- 
grounds by a redefinition of the closed string operators, explicitly at the level of the Virasoro constraints, was 
rediscovered recently [30] . 

^ We suspect that there is a relation with the work of ref. [34] , which recently discussed the organization of 
a number of topological string theory structures in terms of modes of a chiral boson Interestingly, in that 
work, of which we do not understand much (due to no fault of the authors) an exponentiation of the boson 
also plays a natural role in the insertion of D-branes in the background, which reminds us of the observation 
of ref. [20], discussed above. Clearly this connection deserves to be better understood. 
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scaled eigenvalues, and are known from explicit computations (at least for the unitary members 
of the (2, 2k — 1) series [12]) to match the tensions of the D-brane of the continuum theory. The 
interpretation in terms of D-brane world-volume theory and tachyon potentials is discussed 
in detail there. Sections El and d review and complete the derivation of aspects of the 
doubled scaled matrix model origins of the type OA string equations, explains the connection 
to the bosonic models, and discusses the differences in non-perturbative effects. This includes 
a reminder of the role of the two seas of eigenvalues in stabilizing the non-perturbative physics. 
Section |H1 turns to backgrounds with non-trivial open string or R-R flux. The central object 
in this section is the resolvent, which plays two connected roles here. The first is as a probe 
of the matrix model in the language of macroscopic loops, or as the effective force on a scaled 
eigenvalue. We re-derive the effective (tachyon) potential of section |3] from this point of view, 
now including terms higher order in perturbation theory. The second role is as a source in the 
string equations which inserts a number of D-branes, or units of flux, into the background. 
This section has both of these roles working in tandem, since we study the loop operator 
(eigenvalue dynamics) in the presence of the non-trivial backgrounds. So after a brief review of 
the properties of the string equations for these situations, we remind the reader of the properties 
of the scaled loop operator in these theories, and explain the systematics of developing the 
results for the operator perturbatively, using the string equations together with the Gel'fand- 
Dikii equation. We derive compact expressions for the loop operator and hence the eigenvalue 
distributions and effective potentials to disc order (and to annulus and torus order), and analyze 
the results. In section El we remark upon the fluidity of the formalism in its treatment of the 
closed string and open string sectors. We see that the open string (or flux) backgrounds can 
be implemented in terms of a specific preparation of closed string operators. This is made 
much more transparent — and suggestive — in an associated conformal field theory language, 
which organizes everything. We remind the reader of the formalism [26,27] of the Z2 twisted 
boson, and of how the entire content of the string theory is expressed as Virasoro constraints on 
the partition function. This partition function is realized in the twisted boson framework as a 
coherent state. We then discuss how [20] the open string or flux background is simply realized 
by placing a vertex operator into the coherent state background. This shows the remarkable 
transparency of the formalism with regards working with open or closed string backgrounds. 
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2 The String Equations 



Consider the following equations: 

n = 2uTR{z,a) , (1) 

and 

{u - a)n^ - -nn" + ]in'f = . (2) 

2 4 

Here u{z) is a real function of the real variable z, a prime denotes ud/dz, where z/ is a constant, 
and a and T are parameters we shall discuss further shortly. We define the quantity: 



fe=o ^ ^ 



(3) 



where the Rk {k = 0, . . .) are polynomials in u{z) and its 2;-derivatives. They are related by a 
recursion relation: 



4+1 = - uR'k - T^u'Rk , (4) 



^--K K 2' 

and are fixed by the constant Rq, and the requirement that the rest vanish for vanishing u. 
The first few are: 

i?o = ^; Ri^-\u; R2^^{3u^-u") . (5) 

The kih model is chosen by setting all the other ts to zero except to = z, and tk, the latter 
being fixed to a numerical value such that TZ — Vk — z. The are normalized such that the 
coefficient of u'^ is unity, e.g.: 



Vi^u, V2^u --U , Vs^u -uu --^y'^) ■ (6) 



The quantity 



1 92 
i?(^,A) =< ^l-^^^— > , where = -i/^— + 'u(^) , (7) 



is the diagonal of the resolvent of the Hamiltonian H and it satisfies the Gel'fand-Dikii equation 
[35-37]: 

4iu-X)R:' -2RR" + {R'f ^1 . (8) 



The function u defines the partition function Z — exp(— F) of a string theory via: 

(9) 



IJ,=Z 
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where /i is the coefficient of the lowest dimension operator in the world-sheet theory. It is the 
world-sheet cosmological constant in the case A; = 2 for the bosonic string (in which case we 
take equation ([T))) and /c = 1 for the type OA string, in which case we take equation (j2))). For 
the case k = 1 for the bosonic string, /i couples to a boundary operator measuring lengths. 

The equations ((T)) and © have exactly the same physics for large z (equivalent to /i; we will 
sometimes use the two interchangeably), which corresponds to a perturbative regime. Let us 
set r = and cr = for now. For example, for k = 2, we have for the free energy (discarding 
non-universal terms): 

4 1 7 

F=— 0"^ + — logu + (10) 

15^' 24 1440^" ' ^ ^ 

where the dimensionless parameter gg = z//i~^/^ is the closed string coupling. We see the 

expansion in terms of closed string worldsheets of topology x = 2 — 2/;, (/i is the number of 

handles), each coming with factor g~^. The dependence on the sphere displays bosonic KPZ 

scaling F{fi) ~ /i^"'''^*'' for the critical exponent (string susceptibility) 7str = —^/k, and k = 2 

here. For = 1, Z ~ /i'^, which should be discarded with other non-universal physics, and so 

this theory has no contribution from any closed string world-sheets. It is a topological theory, 

as we will recall below. The closed string coupling for this k = 2 case is gs = ufj.~^^'^. For 

general k, it is: 

9s = , (11) 

u{z) = zT^ + --- (12) 

F^^- (13) 

Non-perturbatively, the situations for each equation, and hence the models they describe, are 
different. For the case of equation (^, there is a well-known instability [1-3]. The physics 
is non-perturbatively problematic for even A;, there being no sensible real solutions for u(^z^ 
which match onto the large z perturbation theory just discussed. This fits with the fact that 
the model is rendered unstable through a decay channel mediated by instanton effects [25, 38] 
which we now think of as the nucleation of D-branes. 

Let us turn to equation (0), starting again with a and F set to zero. This equation enjoys 
better non-perturbative properties than the original string equation ((T)), (see refs. [19,39-43]) 
while sharing the same perturbative physics for large positive 2;, as follows from the fact that 
7^ = is also a solution of this equation, and that for large \z\, the leading order equation 



which fits with the fact that 



and hence: 
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is u(u^ — z) = 0. So the behaviour u = -\- . . . can be chosen for large positive z, and the 
perturbative expansion is exactly the same as that given in equation (fTIH) . 

For large negative z, however, we can choose a completely different behaviour for solutions of 
equation Q from that available for equation (P). We can choose the asymptotic m = in that 
regime, having the physics on the sphere vanish for all k. The next non-vanishing order is of 
a universal form for all A;, and is: 

u{z) = --- + ... , (14) 

yielding a universal torus contribution for all models. The full solutions of the equation with 
these asymptotics for large ±2;, for all k, are believed to be unique and pole-free supplying a 
fully consistent non-perturbative physical string theory. See for example the bottom curves of 
the families depicted in figure Elfaj and (h). 




Figure 2: Numerical solutions to equation ^ for u{z): (a) The case A; = 1, and (h) the case k = 2, 
for a range of values of the parameter T, with a = 0. The very bottom curve in each case is the case 
of r = 0. The solutions all asymptote to u = z^^^ for large positive z and u = for large negative z. 

The cases of A; = 1,2 and 3 have been studied extensively, and (at least) A; = 1 is known to 
exist uniquely (since it maps [44,45] to the case of the simplest Unitary matrix model critical 
point [46,47] which defines unique physics [48,49]) and this at least implies a unique family for 
all k connected by KdV flows [42]. 
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These solutions were proposed as non-perturbative completions of the sick bosonic string mod- 
els a while ago [19,39-43]. That statement has now been refined by recent work [21]: The 
models have an interpretation as type OA models. This is consistent with a number of facts, 
among which is the fact [21] that the spectrum of operators of the Liouville dressed (2, 2k — 1) 
conformal minimal models match exactly those of the super-Liouville dressed (2, 4k) supercon- 
formal minimal models. 



3 Matrix Models and Double Scaling 1: Bosonic Strings 

Both sets of models can be derived from one-matrix models in the double scaling limit and 
we can choose Hermitian matrices if we so wish. Equation for F = arose for Hermitian 
matrix models, which can be written in the form: 

Z = y^Mexp |-^Tr (\/(M))| , (15) 

where M is an NxN Hermitian matrix, and the coefficients of the polynomial potential V{M) 
and the parameter 7 are tuned to critical values as 00. Recall that the Feynmann 

diagrams of the matrix model are interpreted as dual to discretisations of the string worldsheets, 
organized in a genus expansion for large A^. The path integral is then to be thought of as a 
particular regularization which enumerates the sum over surfaces. The typical size of a cell in 
the discretization is set by a parameter we shall call 6. The critical potentials are those which 
pick out surfaces with diverging area, so that in the continuum limit 6 = these give finite 
area contributions to the path integral. 

The typical method of study [50] is to diagonalize M using a unitary matrix U into the form 
M = UAW, where A = diag {Ai, A2, . . . , A^v}. There is a Jacobian for this change of variables, 
which is the square of the Van der Monde determinant nj<i('^« "~ -^i)' resulting model 

is (after dividing by a normalization factor Vol{U{N))): 

Z = |j3rfA.n(A.-A,)^exp|-^^V(A.)| , (16) 

An efficient way of working [51] is to use polynomials Pn(A) ~ A" + ■ ■ ■ which are orthogonal 
with respect to the measure (i/i(A) = (iAe~^/'^^^'^\ 

P„(A)P„(A)rf^(A) = hJnm . (17) 
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The polynomial AP„(A) must be equal to P„+i plus a correction, which for even potentials (on 
which we will focus here) is simply: 

AP„,(A) = P„+i(A) + P„P„_i(A) , (18) 

where we have defined the coefficient Rn. Using the orthogonality condition it is given in terms 
of the hnS as Rn+l = hn+i/hn- 

After a little more work, one can write the vacuum energy of the model as [51]: 

N 

F = -Y,{N-n)\ogRn + --- , (19) 

n=l 

where the ellipsis denote parts which are non-universal, and the entire problem is one of de- 
termining the P„s. 

The content of the matrix model can be extracted from the relation: 

/d 
d^l{X)Pn-l{X)—Pn{X) = nhn , (20) 

which amounts, after integrating by parts, to a compact master equation: 

^<n-l\V'{X)\n>='^ , (21) 

where we have used a bra-ket representation of the orthogonal polynomials and the integration 
over product pairs of them with respect to the measure. Functions integrated with respect to 
the measure become operator valued in this representation: 

\n >= , n\n >= n , A\n >= \n — I > , A^\n >= \n + I > 



n 



< m\n >= 5mn , ARn = Rn+lA , A^ Rf, = Rfi-lA^ ■ (22) 

The recursion relation (fTHj) define 

A = ^/R^A^ + A^/R^ , (23) 

which allows a rewriting of the master equation in the equivalent form: 

< n\XV'{X)\n >= j^{2n + 1) , (24) 

which will be useful later on. It is convenient to write our master equation ()21j] a touch more 
compactly as: 

<n|^V(0)|n>=^ , (25) 
11 



where 

= + ^i?- , = Xy/h^ . (26) 

The master equation is now actually a recursion relation for the as can be seen by expanding 
for a given potential V{\). The higher the order of V{X), the more terms there are in the 
recursion relation. 

The continuum limit is then easily approached in these variables. At large N we pass to a 
continuum limit in eigenvalues space as well, and write x = n/N, R{x) = R{n/N), e = 1/N. 
The master equation is: 

<x\A^V'{(i))\x>='yx , (27) 

and the free energy is: 

F = -N'^ I dx{l - x) \ogR{x) . (28) 
Jo 

Focusing on the sphere contribution means dropping terms subleading in N and so the difference 
equations become polynomial relations in R{x): 

1 - 7x = P{R{x)) . (29) 

The required divergent behaviour of the free energy F ~ (7c — occurs when the polyno- 

mial in R{x) acquires k multiple roots: P{R) — {R — R^^. This amounts to choosing particular 
values of the couplings in the potential [52] and also tuning 7 to a critical value 7c. The critical 
value of R{x) at this point can be chosen to be equal to 2, and in the neighbourhood ^x — 1, 
without loss of generahty. We can choose a potential of degree 2k to get the above behaviour, 
and a convenient basis is [53] : 

VM = 

Finite physics is extracted by allowing R{x) and 7a; to approach their critical values at the 
right rate as 5 — > 0. The quantity e = 1/N must also be scaled with 6 so as to approach zero 
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at the correct rate: 



7 = l-^5^^ 

R{x) = 2- u{z)6'^ ; 

(Numerical prefactors of the scaled quantities ^,z,u{z) and u are chosen for convenience.) 
To go beyond the sphere one simply keeps the contributions in the master difference 
equation, and in the large N limit it becomes a differential equation. Inserting the scalings 
into the master equation and expanding about 6 = gives for the first non-vanishing terms, 
a differential equation of order 2/c — 2 at order 6"^^. So in taking the limit 6^0, this is the 
surviving physics and it is in fact equation dH), with F = 0, cr = and all ts set to zero except 
to z and t^- 

By introducing dimensionful coefficients tk for each potential, we can consider the general 
model. From the point of view of the kth theory, the other tkS represent coupling to closed 
string operators Ok- It is well known that the insertion of each operator can be expressed in 
terms of the KdV flows [24, 54] : 

I - • («^) 

The operator Oq couples to to, which is in fact — 4z, the cosmological constant (in the unitary 
model). So Oq is often referred to as the puncture operator, which yields the area of a surface 
by fixing a point which is then integrated over in the path integral. The function u itself can 
be thought of as the two point function of the puncture operator. 

4 Effective Tachyon Potentials and D— Branes: Bosonic 
Strings 

It is also important to keep in mind the eigenvalue picture. The eigenvalues take values on the 
real line. The Van der Monde determinant acts as a repulsive potential, driving them apart, 
but they are confined by the potential V{X) to form a droplet of finite size at large A^. At 
large N, at the spherical level, one can solve for a density p(A) of eigenvalues on the line. A 
relation between the orthogonal polynomial quantity R{x) and the eigenvalue density p(A) can 
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be derived at sphere level [51]: 

dx e{AR{x)-x^) 

' Jo v/4i?(x) - A2 ■ ^ ^ 

From this expression, it is clear that the ends of the eigenvalue density are located at Ac = 
±2a/R^. Recalling that the neighbourhood of a; = 1 is where the interesting physics comes 
from, we must pick one of them to which to scale our physics: 

X = 2V2~^6\ (34) 
v2 

This, together with the scalings gives an expression for the scaled density: 

^^''~~1^J ^A, - uiz) ~ l67 y^v. ^/XT^^ ' ^^^^ 

where we have used that u{z) = z^^'^ on the sphere. Explicit computation of the integral 
yields [40]: 



k\6 



2k-l 



= ^ g ____(,. _ ,t)»-S,.-* + . . . , (36) 

where we remind the reader again that this is valid only at the sphere level. 

Crucially, note that since there is a common factor (A^ — /i^/'^)^, the eigenvalue density vanishes 
at Xs = fi^^'^ as A^. There are k — 1 more zeros which are located away from this endpoint. 
For k even there is another zero at A < and k — 2 zeros off the real line in pairs. For k odd, 
there are no other zeros on the real line, and only pairs of zeros off the real line. In general, 
p(A) is real on the part of the real line where the eigenvalue are located but its complex parts 
are meaningful also, as we shall see. 

A most important quantity is the effective potential that an individual eigenvalue sees [25]. 
It is minus the change in energy as the eigenvalue moves from Aj to Xf, and is given by the 
expression: 

V^sW = N r dX{V'{X) - 2F(A)) = 27rArim f ' dXp{X) , (37) 



where 

F(A)=/#| (3 



is the matrix model resolvent. So starting with the eigenvalue on the cut and moving it to 
anywhere else on the cut costs no energy, since p(A) is real there. In scaled variables the 
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effective potential becomes [40]: 

KSHAs) = r S'piX^K , A. < /ii . (39) 

Inserting our previous expression for the density gives [40]: 

Looking at the first two cases explicitly is illustrative: 

visits) = ;^(/i-A.)i + -- - , A. </i, 

vS\Xs) = |;(/i^-A.)M-3^(/x^-A,)t + --- , K<^^. (41) 
Their behaviours are sketched in figure El 

The key point is that there is a turning point of the effective potential at every zero of the 
eigenvalue density. In particular we see that the even k cases have a maximum of the potential 
at one point on the real axis at A^ < 0. The odd k cases have an effective potential which 
forever rises to the left. To the right in figure El is the sea of the other eigenvalues, starting at 
Xs = jJ^^ and going off to infinity. 

In the spirit of ref. [4], the interpretation of this physics in the continuum string theory is 
as follows. The matrix model, once double-scaled, is the theory on the world-volume of 
unstable D-branes in the continuum string theory. This model is constituted by a matrix- 
valued tachyonic field made by open strings interconnecting the N D-branes. The computation 
above is a precise determination of the potential of the tachyon field at sphere level in string 
perturbation theory. We have a precise model in which we see holography at work, together 
with a model of tachyon condensation. To the right is the sea of eigenvalues. The collective 
dynamics of these represents closed string physics: Bosonic non-critical string theory given 
by Liouville theory coupled to the {2k — 1,2) conformal minimal models. For this model, 
the full family of critical potentials given in equation (jHUj) forms a basis for the closed string 
deformations. They can be added to the model with coefficients and are dual to inserting 
closed string operators Ok to the model. These are the Liouville dressed conformal field theory 
operators. As is well known [24, 54] , this family of operators is organized by the KdV flows of 
equation 

There is more, however [4]. A single eigenvalue can be removed from the sea and moved up the 
potential. This describes the creation of a D-brane of the continuum string theory, and it can 
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Figure 3: The effective potential for one eigenvalue in the scaling limit, computed at string theory 
tree level. This is equivalent to the tachyon potential on the world-volume theory of the D-branes. 
The sea of eigenvalues is to the right in each sketch, on which V^s = 0. The collective physics of the 
sea describes closed string physics. A D-brane is described by a single eigenvalue separated from the 
sea, and its decay by rolling back of the eigenvalue into the sea. The energy released in decay sets 
the tension of the D-brane. Shown is the difference between the D-brane tension in the bosonic (rb) 
and that of the type OA theory (tqa) for the cases oi (a) k odd and (h) k even. In the OA case, the 
potential stops at Ag = 0, setting a lower value for the D-brane tension. 

decay back to closed strings by rolling back down the potential to rejoin the sea. Consistent 
with this identification is the observation of the behaviour of the tension of the D-brane. The 
tension of the D-brane should be the energy needed to move up the potential to a definite 
position which goes like: ~ c/i^, where c is a constant. Upon examination of the potential 
V^^ it can be seen to behave as: 

r^V^^{c^^^)r.C^ = ^, (42) 

^ 9s 

where C is a constant. This is the correct behaviour for a D-brane's tension. 

For k even, there is a special class of D-branes formed by placing an eigenvalue at the maximum 
of the potential, located at = Ac. This is a D-brane which dominates the non-perturbative 
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contributions to the closed string physics, an observation going back to ref . [55] . Its tension is: 



ri"^ v}s\Xc) , for k even. (43) 
For example, the case k = 2, since the potential has a maximum at Ac = — we get [55,56]: 

(2) _ (2) /_1 A ^ _ .44. 

This actually shows up as the action of the D-instanton effects arising by studying the leading 
contribution to effects invisible in perturbation theory. Taking a solution Uq{z) of the equation: 

-^^0 + ul = z , (45) 

one can consider a perturbation e{z), and study another solution u{z) = UQ{z)+e{z). Assuming 
that e is small, and of the form e = exp(— /), where / >> /' >> /", etc., the result [55] is for 
uo = z^': 

This D-brane, (and its analogue for all k even) as can be seen from the effective potential, 
destabilizes the theory for k even, since the eigenvalue can roll off to infinity to the left of the 
figure ini(^&/ This process dooms the k even models [25,38]. Note that we should be careful not 
to confuse the instanton effects derived from analyzing perturbation theory and those of the 
real destabilising instanton see in the effective potential. They match in this case, but existence 
of one does not imply the other, as we shall recall for the type OA models [40]. 



5 Matrix Models and Double Scaling II: Type OA Strings 

The equation (j2I) also arises from matrix models in a very natural manner. It first appeared in 
the context of complex matrix models [39], but can be derived from Hermitian matrix models 
as well. The key point is that the complex matrix models can be formulated as usual in terms of 
eigenvalues. A, of the matrix M^M, which are positive. Therefore the problem naturally places 
a boundary or "wall" at A = 0. The large positive z regime corresponds to the distribution 
pulling way from the wall, and so the system resembles the eigenvalue distributions of the one- 
cut Hermitian matrix models and yields the same perturbation theory. The large negative z 
limit is very different however, as the distribution pushes up against the wall, and gives very 
different physics which we will discuss at length shortly. Another way to look at it is in the space 
A where A^ = A. This is the natural eigenvalue space of one matrix. Then we see that there 
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are two identical eigenvalue distributions, one on A G [0, +00) and the other (a mirror image) 
on A G (—00,0]. See figure IHf^a/ These distributions either pull apart {z +00) or push 
together z ^ —00. For completeness, it is worth recalling how to derive the string equation. 



f ^ff 



(a) 




1^' 



(b) 



Figure 4: (a) Sketch of the effective potential on the covering space for arbitrary k. (b) As fi decreases, 
the eigenvalues fill up closer to the top of the potential, decreasing the tension of the available D-brane, 
until at ^ = 0, there is no brane available. 

We choose the case of an Hermitian matrix model, and we simply introduce a boundary in the 
eigenvalue space. One can define polynomials orthogonal to the measure djjL^X) = dXe~^^"'^^'^^ 
on the interval [— S, E], modifying equation (fT7j) : 



Pn{\)Pra{X)dfi{\) = K^nm ■ 



(47) 



The non-trivial (for even potential V"(A)) relations generalizing equation ()21|1 (or more specifi- 
cally, the form given in equation (jUJ) are [57]: 



<n-l\V'(X)\n > 
< n\XV'{X)\n > 



7?2 27P„(S)P„_i(S) _Ny.j.s 
h -TT , — e T 



7 /,o^ , ,^ ,^Pn(S)2^_^y(s) 



(48) 



18 



where we have used that, for even potentials, P„(— A) = (— 1)"P„(A). We can ehminate the 
boundary terms to give: 

' < n\\V {X)\n > ) ( — <n-l\\V {\)\n-l> 

< n - l\V' {\)\n >\ . (49) 



4 V ' ' / V 

7?T, 



njr: 



Note that 



< n\XV {X)\n >= ^Rn < n\A^V {X)\n > +^Rn+i < n\AV {\)\n > , (50) 
and move to the continuum variables as before, defining 

il{x) = -fx- < X - e\V'{(i))\x > . (51) 

Then the equation becomes: 

[n{x) + n{x + e)] [n{x - e) + n{x)] = ^^^^^i^) ■ (^2) 

Using the scalings given in equation (jSH), the critical potentials given in equation (jHUjl . and 
placing the walls at the ends of the critical eigenvalue density at S = _Rc we get the string 
equations (j21) in the limit 6^0, having appeared at first non- vanishing order 5'^'^+^. 

6 Matrix Models and Double Scaling III: Including Open 
Strings 

In this section, we discuss how the full string equations and (j2|), with all parameters switched 
on, arise in the matrix models. 

Let's first turn to non-zero a appearing in equation Q. In fact, a naturally appears in the 
problem as the scaled position of the wall, which should be a physical parameter as well [57]. 
In taking the double scaling limit, scaling the position according to S = Rc{l — a6'^/2), yields 
the equation ^ with non-zero a. Well behaved solutions for non-zero sigma are known to 
exist, and their behaviour is given numerically in figure El Next, let us note that there is an 
important quantity which arises in the continuum limit of the matrix models which will play 
an important role here and later. In the bra-ket formalism introduced earlier, the operator A 
in equation ()2(i|l can be directly formulated as a differential operator [53]. This comes because 
the shift operators ^, get represented as derivative operators exp (ec^a;) , exp (— e^^:) acting 
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Figure 5: Solutions to equation © for u{z), with asymptitics given in figured (a) The case k = 1, 
and (b) the case k = 2, for a range of values of the parameter a with F = 0. The middle curve in each 
case is the case of o" = 0. 



on functions of x = n/N. Recall also that there is a specific scaling for the objects from which 
A is made, given in equation (jHT|) . so we can expand it and see what its scaling is in the limit 
6—^0 and we find that: 

A = + ^ (u^-^ - u{z)] 52 + . . . = A, - + ■■■ (53) 



V2 

The Schrodinger operator Ti. plays a crucial and very natural role in what is to come. 

The matrix model origin of the string equation ([Q) with non-zero F arises by adding a term 
[17, 18] 

AV{M) = ^Tr log (1 - s^M^) (54) 

to the matrix model potential. In the world-sheet dual to the Feynmann diagrams, this cor- 
responds to adding holes of all sizes, by expansion of the logarithm. This new term adds new 
pieces to the master equation (j2H): 



V l-s2A2 

On the sphere, there is still the critical behaviour where large area surfaces dominate, near 
R = Rc, but there is an additional critical point near a critical value of s set by = 
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2a/R^. This arises because, e.g., on the sphere, the last term in equation reduces to 
-{T/N){1 - 4s2i?(x))-i/2. So the lengths of loops can chosen to diverge if we send s — > Sc, but 
give finite length loops as 5 —> if we scale at the right rate. Since near the critical point the 
new terms are: 

- — r , (56) 

N[A{.sl-.s^)R, + iR,-R)/R,]'2 

this rate turns out to be = — a6^. The factor of scaling as ^^'^+1 produces the 
appropriate damping of the divergence and the string equation ((H) results (as before) at order 
S'^'^, the resolvent arising from the appearance of Ti. in the scaling part of A, appearing at the 
same order in 5 as a does. 

The same procedure can be applied to the master equations ()49|1 arising from placing walls 
in the system. As happens before, order ^^'=+2 is where the non-trivial physics appears, and 
the lefthand side of the string equation Q appears at this order, while the right hand side, 
z/^r^ appears multiplying —1 times the left hand side of the Gel'fand-Dikii equation (jH)), with 
A set to a. Thus we arrive at equation (j2)). So this is how the position of the walls and the 
weight factor of the loops become related in the continuum model. In effect, the argument of 
the logarithm diverges (so as to give finite length loops in the continuum limit) at the same 
place in eigenvalue space at which the walls are located. 

Another way to see that this must be true is to imagine that equation Q was true, but for 
another parameter r, instead of a. How can we fix r? Taking sums and differences of the 
equations we see that we have a condition that either t = a and the equation Q follows, or 
we have 7?. = 0, in which case F = 0. This matches what we learned from the matrix model. If 
we do not scale the parameters so as to make the logarithm diverge precisely at the ends, then 
it produces no new physics. 



7 Effective Tachyon Potentials and D— Branes in Type 
OA String Theory 

Turning to the issue of instantons arising from analyzing perturbation expansions, we note [40] 
that there is another instanton contribution to the large positive z perturbation theory. Search- 
ing as before for the perturbation e = e~'^, we get the same instanton as for the equation (^, 
and in addition, there is an equation: 

K-^(/')^)7^^ = 0, (57) 
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which, upon using Uq = z^^ gives instanton: 

4fc 1 

e ~ e 9^ . (58) 

This instanton can be connected to the physics of the effective potentials we have already 
studied. See figures El and IHf^a/ For large positive z, the results from before will still apply 
(after a shift of the origin of coordinates to place one of the walls there), modified only by the 
fact that the eigenvalues must be positive. Looking at the effective potentials, we see that since 
there is a wall at = 0, the instanton at negative previously encountered for k even cannot 
destabilize the system. There is an instanton for all k at the wall, and its tension is given by 
V^g''(0), which a computation shows is: 

4-k ii^^'ik 

.„, = l.,(A, = 0) = — ^i^. (59) 

This is the same action as the new instanton perturbations of the solutions we saw above. So 
we see that although the same instanton which destabilised the bosonic theory is present in the 
analysis, it does not represent an instability [40]. 

Since the role of a in the matrix model is to move the wall, as we saw in the previous section, 
it can be seen as an interpolating parameter in some sense, since the limit of large negative a 
approaches (but never reaches) the physics of the unstable bosonic models [58]. See figure 
Notice for example how as a becomes more negative, for the case k = 2, the system begins 
to develop more fluctuations at negative z, ultimately heralding the poles that exist on the 
negative z axis for the Painleve I equation (equation ([[J for k = 2.). Studies of this sort were 
presented in ref. [58]. 

The physical role of a from the point of the continuum string theory is more subtle. It represents 
a non-perturbative contribution to the coupling of a boundary operator Ob which measures 
the length of world-sheet boundaries or loops [58]. It is a boundary cosmological constant and 
we will examine it in some more detail in the next section. Note that for the models defined 
by equation with F = 0, the boundary length operator for the kth model is identified with 
Ok-i- See ref. [59]. In particular, for = 1, the operator is Oq which couples to to ~ z. For 
this case, (and for the k = 1 solution of equation Q for large positive z) it is interesting to 
note that a can be absorbed into z without changing the universal physics. So we see that a 
is indeed a natural realization of the boundary cosmological constant, by matching onto the 
perturbatively identified boundary operator. 

Recall that the k = 1 model has no expansion in terms of closed string world-sheets. The the- 
ory's entire content comes alive when the other operators Ok are inserted, and their correlation 
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functions at genus h are then seen to be isomorphic to intersection numbers on the moduh space 
of Riemann surfaces of genus h: it is a topological theory. The k = 1 solution of equation Q 
(for large z) shares this behavior, but develops new features in the non-perturbative regime. 
In fact, this is a precise example of a non-perturbative completion of a topological theory, and 
deserves further study from this perspective. It can be rewritten as the type OB model in the 
same background, as identified in ref. [21] using the change of variables of refs. [44,45]. 

As we will recall below the partition function of the k = 1 theory (exactly for equation (0) 
and at large positive z for equation Q) only has a topological expansion in the presence of 
D-branes, giving only world-sheets with boundaries in the expansion. 



8 Loops and D— Branes for Bosonic and Type OA Strings 

Now we turn to the physics of non-zero F. A study of equation (j21) for non-zero F shows 
that large positive z (or /i) perturbation theory has an interpretation as including world-sheet 
boundaries, and in fact F counts the number of D-branes in the background, as can be seen by 
the fact that a factor F* appears for a surface with the number of boundaries equal to b. For 
example, for k = 2: 

= T.^Kkg's'-'^'rK (60) 

h,b 

The equation supplies fully non-perturbative completions of this physics, including the following 
behaviour at large negative z (or fi): 

(4r2 - 1^ 



4 



log/x (61) 



+ (4F2 - 1)(4F^ - 9)(4F2 - 25) l^^gt + ^(48F^ " 1240F2 + 8371)^^ | + ■ ■ ■ , 

where the first term is universal to all k. At face value, the interpretation in terms of boundaries 
can be made in this regime as well, but explanation for the number of exactly vanishing orders 
in perturbation theory (all surfaces with an odd number of boundaries, for example) is needed. 
Instead, ref. [21] pointed out that in this regime, there are no branes, but instead R-R fluxes, 
using facts from the continuum type theories and the matrix model. Here, we'll study the 
matrix models further and uncover the effective potential beyond tree level in both regimes. 
Its behaviour in the different regimes clearly give support to this interpretation. 
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The large positive z perturbation theory of the equation matches perturbation theory obtained 
from equation (0), which for non-zero F and a was formulated in refs. [17,18]. F appears 
multiplying the resolvent, which in turn satisfies equation (jHJ. In fact, the most efficient way 
to construct a perturbative expansion for u(^z) is to treat equation as an equation for R and 
then substitute it into the resolvent equation, yielding an equation for 71. and hence for ui^z). 
This of course yields the equation (j2I). Of course, this larger equation allows solutions for u{z) 
which are not solutions of equation (^J). These are the solutions which have the well defined 
large negative z asymptotics and accompanying pleasant non-perturbative behaviour which we 
now associate to type OA. See the families of curves displayed in figure El 

That the resolvent appears explicitly in the description of D-branes in the backgrounds is 
interesting because the same object is known to control the behaviour of macroscopic loops in 
the continuum theory [24]. This is worth recalling in some detail: In the matrix model, it is 
natural to consider the Wilson loop: 

W[i) = (^Tre^^) . (62) 

Its Laplace transform is 

1 1 

^^^^ = ^N^'X^^ = I ^^^■'''^W • (63) 
and in the continuum limit it is: 

W{\) = [ . (64) 



The saddle point equation of the matrix model is: 



where P denotes the principal part and the integral is over the support C of the density function, 
p(A), of eigenvalues A. The density is normalized according to J dXp{X) = 1. The real part of 
W is equal to V /2 on the support of p(A) where its imaginary part is non-zero. This is nicely 
expressed by the equation (see ref. [60] for a useful presentation): 

Im[P(A)] = lm{V'W - W^) = {V - 2ReW)lmW = (66) 

with solution 

iy(A) = ^±^vW^4P(A) . (67) 
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where -P(A) is fixed by the asymptotic W{X) ^ — 1/A as A ^ oo, given the normahzation of 
the density. By definition, p(A) is given by the discontinuity of W{X) on the cut C. 

The connection with the continuum loop operator and the continuum resolvent operator comes 
from the definition: 

WiX) = ^ + 6''''w{Xs) (68) 
where, as before, A^ is the scaling part of A. 
Now the marvellous thing [24] is that 



w{Xs) = ldzR{z,Xs)= I dz < z\j^p^—^\z > , (69) 



and it is the Laplace transform of the continuum loop 



oo 



w{e) = J dz < z\e~'^\z > ; w{Xs) = J dee-^''w{e) (70) 

The Schrodinger operator 7i arose earlier in equation ()53|) . Loops correspond to expectation 
values of traces of powers of matrices TrM^ where P = |- This corresponds to cutting a p-sided 
hole in the dual picture of discretized worldsheets. We hold i fixed by sending p to infinity as 
we send 6 —>■ 0. In the bra-ket formalism this amounts to computing the Laplace transform of 
the resolvent R{z, Xg), and we integrate over z (up to some reference value, /i although we will 
still use z for a lot of what follows) to get the expectation value w(£) ~ J dz < z\X^z >, which 
gives the first expression in equations (fTOj) . 

Notice that w{Xs) is also (up to a factor z//2) the first As-derivative of the effective potential 
Ves{Xs), and this gives an efficient way of deriving the effective potentials for the system to 
arbitrary order in perturbation theory, by use of the string equations in combination with the 
resolvent equation. In particular, the procedure is to solve either string equation for u{z) to 
desired order in perturbation theory, keeping some non-zero value of a. This is now to be 
treated as a potential for the Hamiltonian Ti, and one can solve for its resolvent R{z, Xg) to 
that order in perturbation theory by using equation (jH}. The next step is to integrate once 
with respect to z to get the loop operator w{Xs), and one can integrate once with respect to A^ 
to get the effective potential Ves- 

^-^ = IJ wi^s)dXs. (71) 

In fact, it is more interesting to examine the properties of w{Xs), for its zeros tell us about 
the possible location of D-branes. We already know the structure of the zeros at the level of 
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the sphere in perturbation theory, where T is invisible. This procedure allows us to efficiently 
examine the effect of non-zero F on the structure of the zeros. It also allows us to study the 
large negative z behaviour of the effective potential rather clearly. The previous methods were 
limited to the sphere, and in particular the leading behaviour on the sphere vanishes for large 
negative z. 

To start with, take the case k = 1. For large positive z, we have, to disc level: 

u{z) = z+ / + ■ ■ ■ (72) 
- cr 

This can be substituted into equation (jS)), and solve for the resolvent to the same order in 
perturbation theory. The result is: 

R(z, AJ = r ^ r + ■ ■ • , 73 

^ ^ 2(,_A,)i Hz-X,)i{z-a)h ^ ' 

and this results in, after integrating with respect to z: 



w{\s) = (/. - A.)^ + \ \ +■■■ , (74) 

^ (A, - (T)(/i - Xs) - 



and an integration with respect to A^: 



Kff(A.) = T^(/x - A,)i + 2Ftanh-i ( ] _ (75^ 



3z/ \ /i — o" 

Look at the zeros of w{X), we find that we must solve the equation: 

(A, -a)(/i-A,) + ^z/F(/i-(T)5 =0 . (76) 

This is a cute equation. At the level of the sphere, (F = or z/ = 0), it is a quadratic, and we 
see that we have the roots A^ = a, fi. This is what we had before, recalling that a marks the 
position of the wall, and fi is where the eigenvalue density p{Xs) ends. Note that the value of 
the effective potential at fi is zero, while at A^ = a it is 

Kff(a) = ^ ' . (77) 

Note that at cr = 0, this gives tqa- Beyond the sphere we see that F deforms the quadratic, 
changing the location of the zeros: 



a + jj, 1 



As = ^ ± 2 V - Z^)' + 2^rv/^ . (78) 
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To get a feeling for the effects of the parameter F, let us work at large fi, in which case w{Xs) 
becomes: 

w{Xs) = - + + -— , (79) 

^ /i2 2 As - a An 



with zeros at approximately: 



1 z/F , , 1 z/F , , 

A;=a--— , A+ = /i+-— . (80) 



The zero to the right of /i lies on a cut which runs to infinity, marking support of the eigenvalue 
distribution. The other zero can be thought of as a cut which has collapsed to zero size, in 
the language of ref. [21]. As explained in ref. [21], there is a collapsed Riemann surface of the 
form y"^ = f{Xs) where y = w{Xs), extending to the complex A^-plane. The vanishing of w{Xs) 
controls the location of the cuts on the complex plane A^, and the surface is the double cover 
of the plane joined through the cuts, giving a surface of specific topology. We can perform the 
integral (I/ttzz/) f w{Xs)dXs along a closed contour running between the cuts on one sheet and 
closing the loop by going along the other sheet. The normalization is chosen to compute the 
change in the effective potential Vcs, following equation ()4()|1 . For our problem in hand, we see 
that the result comes from the pole at A^ = a, and the result is precisely F. This counts the 
number of branes in the problem. 

We've computed a pleasant general structure for the higher k cases. The leading contribution 
(up to disc order) to u{z) is: 



1 z/F 

fc — 1 1 

z~>^(zk — cr)2 



u{z) = z" + j:^^r—i rr + • • • 



and some algebra produces the resolvent: 



11 1 z^r 

■^[Zk — As)2 (Zk — Xj2(Zk —a)2Z k 



Integrating this up gives 



^ '^"-1 - 1 , „ 1 vTizi -a)^ 



and 
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where we reproduce the previous effective potential expressions on the sphere, but have in 
addition a disc level contribution. 

The location of the zeros of w^''\\s) give the extrema of the effective potential. The issue boils 
down to a (/c + l)th order polynomial in A^, which on the sphere has zeros at a and /i^ for 
positive Xs, and for odd k there are no other real ones (only k — 1 imaginary ones). For even k 
there are k — 1 imaginary ones and a single extra one at — c/i^. At the level of the disc these 
zeros get deformed, as before, but notice that there is again a pole at A = a, which will give a 
non-vanishing contribution of precisely F to the B-cycle integral {l/Triu) f w{Xs)dXs between 
the cut to the right of A^ = /i^ and the collapsed cut located just to the left of A^ = a. 

There are presumably means of deforming the problem in order to have non-zero B-cycle 
integrals between the cut and the other {k — 1) zeros, and these would define other branes too. 
It is not clear whether this can be done in this matrix model description. 

We can just as easily study and construct the same quantities in the negative z limit too. In 
fact, we can write everything succinctly for all k to the first non- vanishing order. We have in 
the large negative z direction: 



z/2 (4F2 - 1) 



Putting this into the resolvent equation (jHJ, we get 



From this we get: 



and 



^1 1 _ 1 z/2(4F2 - 1) 

'^~2{a- A,)i ~ 16 (a - A)f - a^')' ^ ' " 



, 1 a 1 z/2(4F2-1) ^ ^ 

wiXs) = - , ^,. +777: rnn —7^ (§7) 



2(a-A)5 lQ{a-X)Hii 



Kff(AJ =/i(a- AJ^ - — . (88) 

The most interesting thing here is that the zeros are at: 

z/2(4F2-l) ^ ^ 

As = cr - -^-^ { . 89 

There is only one zero. On the sphere, it is at a. The effective potential vanishes there. So 
there is no sign of anything that we would interpret as a sector representing D-branes of finite 
tension. This is consistent with the interpretation that there are no D-branes in this limit. All 
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we have is the eigenvalue sea. In a sense it has filled up beyond the top of the potential and so 
there are no branes. This is what is depicted (somewhat schematically) in figure ). 

Beyond the sphere does not change this conclusion. There is a single zero, and the value of 
the potential does not have a scaling normally associated to D-branes. It is important to 
note that this is because of how equation Q and the resolvent equation works perturbatively. 
First note that the resolvent equation (jHJ cannot itself generate odd powers of Qs or F in any 
expansion for R{z,X). It must get it from the potential u{z) which is governed by the string 
equation. It can do this from (for example) the disc term in the large positive z expansion. 
Equation 121 however, gives only even powers of Qs and for F in a large negative z expansion for 
u{z). This then ensures that the resolvent will never have odd T or Qs, and therefore the effective 
potential V^fr(As) will never give rise to energies which have such dependence. Therefore there 
can be no objects which have their mass-energy originating at disc order in perturbation theory, 
and hence no D-branes for large negative z. This constitutes a proof of this directly in this 
language, and lends further credence to ref. [21] 's interpretation of F as controlling closed string 
backgrounds, namely R-R fluxes. 

It is further interesting to note that if we include the next order in perturbation theory beyond 
the sphere, the value of the effective potential at the zero is: 



which scales as g^, which is the correct scaling for R-R flux in closed string perturbation theory, 
or for a torus contribution. 

Finally, it is interesting to study a family of exact solutions to equation Q for each k. The 
solutions are: 



For these, the analysis above goes through as before, with similar conclusions. It is not known 
what the signiflcance of such solutions might be in this context, but their existence is worth 
remarking upon, especially since (as noted in ref. [40]) all the solutions asymptote to them at 
large negative z, and so they may capture some universal aspects of the physics in this regime. 
We make a further observation about these solutions in the next section. 




(90) 




1/2 (4F2 - 1) 




m = . . . k . 



(91) 
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9 Loop Equations and Virasoro Constraints 



As we've seen, the basic organizing object in the theory for much of our considerations has been 
the loop operator w{Xs), which is the first integral of the diagonal of the resolvent R{z, Xg). We 
have studied it in the presence of D-branes (and fluxes) in the background, which enters the 
matrix model string equations in a natural way, again in terms of the resolvent R{z, a), now in 
terms of the fixed background parameter a and multiplied by F. 

It is natural that the loop expectation value w{Xs) is controlled by the same object that we use 
to place D-branes (and fluxes, for type OA) into the the background, and one might wonder 
if this relation can be made more precise. Can one express the open string background as a 
specific combination of the existing closed string operators? If so, this would be a clear example 
of the interchangeability of closed and open string descriptions of a given background, perhaps 
the sharpest demonstration of open-closed duality. In fact, it can be done [20,28], and it may 
be illustrative of the correct language to use for other examples. 

Let us describe this using the equivalent microscopic loop operator language. The object w{£) 
can be expanded in terms of the microscopic scaling operators 



n=0 V ^ / 

and hence the Laplace transformed operator can be written: 

w{x)^J2^~''~''^n- (93) 

n=0 

Insertion of operator On is equivalent to differentiation with respect to the coupling 

<0„>=.^|;. (94) 

It was shown long ago [26, 27] that the purely closed string theory (F = 0, cr = 0) is equivalent 
to the following constraints on r = = ^yexp(—F)■. 

L„T = , n > -1 , (95) 
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where 



k=i ^ ^ 

J,— n ^ / 



^0 = > ; u + ^ 1 + ^ 



fc=0 
oo 



16 ' 



dtk-ldtn-k 



(96) 



which form a Virasoro algebra: 

[L„, Lrr,] = {n- m)Ln+m ■ (97) 
Note that the L_i constraint, after taking a derivative, is equivalent to the string equation 

v^^f, l\ OF' z 



leading to 

°° / 1 \ 

" 4-^ = 7^ = 0, (99) 



where we have used the first integral of the KdV flows ()H2|1 . in the form: 

dF' 



dtk-i 



Rk . (100) 



The structure of these equations become even more natural when one expresses it in conformal 
field theory language. There is a Z2-twisted boson with mode expansion: 

a^(x) = 5^a„,_.x— i . (101) 

n 

The Virasoro generators are then modes of the stress tensor of this boson, acting on a state for 
which the creation and annihilation operators act by multiplication and derivation: 

Such states are coherent states: 

^ (103) 



< t\ =< 0| exp tnCin+^ 
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and the partition function defines a special state fit in the Fock space of (p{x): 

T=<t\nt> . (104) 

It is clear now from equation (jHSI) and equation (jlUlj) that our loop operator is the twisted 
boson [26]: 

w{x) - 0\dyD{x)\Qt > , (105) 
the square root branch cut we studied earlier realizing the Z2-twist. 

The real bonus comes when we see how naturally the non-zero cr and F case fits into this 
formalism [20]. The microscopic loop equations are: 

Ln ■ n^,T} = , (106) 

where ^ 

L„ = L„ - (1 + n)^a" - . (107) 

In particular, if we use that the resolvent R^^z has an expansion in terms of the Gelfand-Dikii 
polynomials as 



oo 



k=0 ^ ^ 

the first constraint is equivalent to equation (Q), (the rest following from the KdV flows) if the 
following relation is true: 

;|-2r.E^--J^. (109) 

n=0 



which is equivalent to [20, 28] : 



tfc = 4 + 7^ n^^^^- (110) 

(^+2) 

This relation is a preparation of a specific combination of closed string operators, shifting the 
background. But from our studies in earlier sections, we know that this background corresponds 
to D-branes in an expansion in large positive /i, and fluxes for type OA in the large negative /i 
limit. 



In fact, the Virasoro constraints p06p with shifted Virasoro operators ()107p are precisely those 
of a vertex operator of weight F^/4, =: e~ v^'^'-'^^ : and the pleasingly simple relation can be 
derived [20]: 

r{,,r}=<t|:e-^'^('^^|fif> . (Ill) 
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Finally, note that the weight of the vertex operator is r^/4. The weight of the twisted boson it- 
self is 1/16. This means that when F = 1/2, the combined wavefunction given in equation pilj) 
comprising the partition function will have vanishing weight, as can be seen by looking at the 
n = case of equation ()107|) and (|^. We do not know the significance of this, but it is of note 
that for this case there is a special exact non-perturbative solution of the string equation (j2)) for 
any value of k, as displayed in equation ()91|) . It would be interesting to explore the role of this 
solution further. Note that, for given k, there are k + 1 solutions of this type, with F = m + 1/2, 
m = 0, . . . , k, and the weight of the partition function then reduces to the interesting quantity, 
m{m + l)/4. 



Acknowledgments 

CVJ wishes to thank James E. Carlisle for many questions and discussions. 



References 

E. Brezin and V. A. Kazakov, "Exactly Solvable Field Theories Of Closed Strings," Phys. 
Lett. B236 (1990) 144-150. 

M. R. Douglas and S. H. Shenker, "Strings In Less Than One-Dimension," Nucl. Phys. 
B335 (1990) 635. 

D. J. Gross and A. A. Migdal, "Nonperturbative Two-Dimensional Quantum Gravity," 
Phys. Rev. Lett. 64 (1990) 127. 

J. McGreevy and H. Verlinde, "Strings from tachyons: The c = 1 matrix reloated," 
|hep-th70 304224. 

V. Fateev, A. B. Zamolodchikov, and A. B. Zamolodchikov, "Boundary Liouville field 
theory. I: Boundary state and boundary two-point function," ,hep-th/0001012, 



J. Teschner, "Remarks on Liouville theory with boundary," hep-th/0009138 



A. B. Zamolodchikov and A. B. Zamolodchikov, "Liouville field theory on a 
pseudosphere," |hep-th/0101152 



A. Sen, "SO (32) spinors of type I and other solitons on brane- antibrane pair," JHEP 09 
(1998) 023,pp-th/980814Tj 



33 



[9] A. Sen, "Rolling tachyon," JHEP 04 (2002) 048, |hep-th/02032T7 

[10] I. R. Klebanov, J. Maldacena, and N. Seiberg, "D-brane decay in two-dimensional string 
theory," JHEP 07 (2003) 045, ,hep-th/0305159t 

[11] J. McGreevy, J. Teschner, and H. Verlinde, "Classical and quantum D-branes in 2D 
string theory," |hep-th/0 305194< 

[12] S. Y. Alexandrov, V. A. Kazakov, and D. Kutasov, "Non-perturbative effects in matrix 
models and D-branes," JHEP 09 (2003) 057, ,hep-th/0306177| 

[13] M. Gutperle and P. Kraus, "D-brane dynamics in the c = 1 matrix model," Phys. Rev. 
D69 (2004) 066005, hep-th/0308047l 



[14] E. J. Martinec, "The annular report on non-critical string theory," hep-th/0305148 



[15] V. Schomerus, "Rolling tachyons from Liouville theory," hep-th/0306026 



[16] A. Sen, "Open-closed duality: Lessons from matrix model," hep-th/ 0308068^ 



[17] V. A. Kazakov, "A Simple Solvable Model Of Quantum Field Theory Of Open Strings," 
Phys. Lett. B237 (1990) 212. 

[18] I. K. Kostov, "Exactly Solvable Field Theory Of D = Closed And Open Strings," Phys. 
Lett. B238 (1990) 181. 

[19] S. Dalley, C. V. Johnson, T. R. Morris, and A. Watterstam, "Unitary matrix models and 
2-D quantum gravity," Mod. Phys. Lett. A7 (1992) 2753-2762, ,hep-th/9206060, 

[20] C. V. Johnson, "On integrable c < 1 open string theory," Nucl. Phys. B414 (1994) 
239-266, |hep- th/930 1112| 

[21] I. R. Klebanov, J. Maldacena, and N. Seiberg, "Unitary and complex matrix models as 
1-d type strings," ,hep-th/0309168, 

[22] C. V. Johnson, "Non-perturbative string equations for type OA," JHEP 03 (2004) 041, 
|hep-th/0311129[ 

[23] C. V. Johnson, T. Morris, and B. Spence, "Stable nonperturbative minimal models 
coupled to 2-D quantum gravity," Nucl. Phys. B384 (1992) 381-410, |hep-th/9203022| 



34 



[24] T. Banks, M. R. Douglas, N. Seiberg, and S. H. Shenker, "Microscopic And Macroscopic 
Loops In Nonperturbative Two- Dimensional Gravity," Phys. Lett. B238 (1990) 279. 

[25] F. David, "Loop Equations And Nonperturbative Effects In Two- Dimensional Quantum 
Gravity," Mod. Phys. Lett. A5 (1990) 1019-1030. 

[26] R. Dijkgraaf, H. Verlinde, and E. Verlinde, "Loop equations and Virasoro constraints in 
nonperturbative 2-D quantum gravity," Nucl. Phys. B348 (1991) 435-456. 

[27] M. Fukuma, H. Kawai, and R. Nakayama, "Continuum Schwinger- Dyson Equations And 
Universal Structures In Two-Dimensional Quantum Gravity," Int. J. Mod. Phys. A6 
(1991) 1385-1406. 

[28] Y. Itoh and Y. Tanii, "Schwinger-Dyson equations of matrix models for open and closed 
strings," Phys. Lett. B289 (1992) 335-341, hep-th/9202080 

[29] M. Fukuma and S. Yahikozawa, "Comments on D-instantons in c < 1 strings," Phys. 
Lett. B460 (1999) 71-78, ;hep-th/9902169| 

[30] D. Gaiotto and L. Rastelli, "A paradigm of open/closed duality: Liouville D-branes and 
the Kontsevich model," |hep-th/0312196] 



[31] N. Seiberg and D. Shih, "Branes, rings and matrix models in minimal (super) string 



theory," JHEP 02 (2004) 021, hep-th/0312170 



[32] D. Kutasov, K. Okuyama, J. Park, N. Seiberg, and D. Shih, "Annulus amplitudes and 



ZZ branes in minimal string theory," hep-th/ 0406030, 



[33] V. A. Kazakov and I. K. Kostov, "Instantons in non-critical strings from the two-matrix 
model," |hep-th/0 403152| 

[34] M. Aganagic, R. Dijkgraaf, A. Klemm, M. Marino, and C. Vafa, "Topological strings and 
integrable hierarchies," |hep-th/ 0312085| 

[35] I. M. Gel'fand and L. A. Dikii, "Fractional Powers of Operators and Hamiltonian 
Systems," Funct. Anal. Appl. 10 (1976) 259. 

[36] I. M. Gel'fand and L. A. Dikii, "The Resolvent and Hamiltonian Systems," Funct. Anal. 
Appl. 11 (1976) 93. 



35 



[37] I. M. Gel'fand and L. A. Dikii, "Asymptotic behavior of the resolvent of Sturm-Liouville 
equations and the algebra of the Korteweg-De Vries equations," Russ. Math. Surveys 30 
(1975) 77-113. 

[38] S. Dalley, "Instability of even m multicritical matrix models of 2-D gravity," Phys. Lett. 
B253 (1991) 292-296. 

[39] S. Dalley, C. V. Johnson, and T. Morris, "Multicritical complex matrix models and 
nonperturbative 2-D quantum gravity," Nucl. Phys. B368 (1992) 625-654. 

[40] S. Dalley, C. V. Johnson, and T. Morris, "Nonperturbative two-dimensional quantum 
gravity," Nucl. Phys. B368 (1992) 655-670. 

[41] S. Dalley, C. V. Johnson, and T. Morris, "Nonperturbative two-dimensional quantum 
gravity, again," Nucl. Phys. Proc. Suppl. 25A (1992) 87-91, hep-th/9108016| 

[42] C. V. Johnson, T. R. Morris, and A. Watterstam, "Global KdV flows and stable 2-D 
quantum gravity," Phys. Lett. B291 (1992) 11-18, hep-th/9205056, 

[43] C. V. Johnson, "Non-Perturbatively Stable Conformal Minimal Models Coupled to Two 
Dimensional Quantum Gravity". PhD thesis, Southampton University (UK), 1992. 

[44] T. R. Morris, "2-D quantum gravity, multicritical matter and complex matrices," . 
FERMILAB-PUB-90-136-T. 

[45] T. R. Morris, "Multicritical matter from complex matrices," Class. Quant. Grav. 9 
(1992) 1873-1881. 

[46] V. Periwal and D. Shevitz, "Unitary Matrix Models As Exactly Solvable String 
Theories," Phys. Rev. Lett. 64 (1990) 1326. 

[47] V. Periwal and D. Shevitz, "Exactly Solvable Unitary Matrix Models: Multicritical 
Potentials And Correlations," Nucl. Phys. B344 (1990) 731-746. 

[48] C. Crnkovic, M. R. Douglas, and G. W. Moore, "Physical solutions for unitary matrix 
models," Nucl. Phys. B360 (1991) 507-523. 

[49] A. Watterstam, "A Solution to the string equation of unitary matrix models," Phys. 
Lett. B263 (1991) 51-58. 



36 



[50] E. Brezin, C. Itzykson, G. Parisi, and J. B. Zuber, "Planar Diagrams," Commun. Math. 
Phys. 59 (1978) 35. 

[51] D. Bessis, C. Itzykson, and J. B. Zuber, "Quantum field theory techniques in graphical 
enumeration," Adv. Appl. Math. 1 (1980) 109-157. 

[52] V. A. Kazakov, "The Appearance Of Matter Fields From Quantum Fluctuations Of 2-D 
Gravity," Mod. Phys. Lett. A4 (1989) 2125. 

[53] D. J. Gross and A. A. Migdal, "A Nonperturbative Treatment Of Two-Dimensional 
Quantum Gravity," Nucl. Phys. B340 (1990) 333-365. 

[54] M. R. Douglas, "Strings In Less Than One-Dimension And The Generalized K-D- V 
Hierarchies," Phys. Lett. B238 (1990) 176. 

[55] S. H. Shenker, "The Strength of nonperturbative effects in string theory,". Presented at 
the Cargese Workshop on Random Surfaces, Quantum Gravity and Strings, Cargese, 
France, May 28 - Jun 1, 1990. 

[56] F. David, "Phases of the large N matrix model and nonperturbative effects in 2-d 
gravity," Nucl. Phys. B348 (1991) 507-524. 

[57] S. Dalley, "On loop equations in KdV exactly solvable string theory," Mod. Phys. Lett. 



A7 (1992) 1263-1272, hep-th/9111064 



[58] C. V. Johnson, T. R. Morris, and P. L. White, "The Boundary cosmological constant in 
stable 2-D quantum gravity," Phys. Lett. B292 (1992) 283-289, hep-th/9206066l 

[59] E. J. Martinec, G. W. Moore, and N. Seiberg, "Boundary operators in 2-D gravity," 
Phys. Lett. B263 (1991) 190-194. 

[60] Y. Makeenko, "Methods of contemporary gauge theory,". Cambridge, UK: Univ. Pr. 
(2002) 417 p. 



37 



